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A simpler and improved utility approximate point scattered function for thin-film converters currently used 
in neutron photographic devices is proposed as a correction method to produce clearer, more realistic images. 
The validity of the model was demonstrated through a simulation experiment. Based on the results, an error 
analysis was carried out, certain corrections were made to the original model, and the final model achieved a 
very low relative error in the simulation experiment. The model can also be optimized for quantitative neutron 
photographic analysis using iterative algorithms to obtain realistic neutron photographic images more quickly. 
At the end of the article, the model is extended to consider the case of energy spectrum hardening by introducing 


a temperature correction parameter. 
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I. INTRODUCTION 


Neutron radiography is a method for surface or internal di- 
agnosis of a sample to be measured using the properties of 
neutron attenuation and scattering. Compared to other ra- 
diographic techniques, neutrons are more complex because 
of the manner in which they interact with the nucleus, and 
the macroscopic cross section does not have a definite re- 
lationship with the nuclide species. These properties pro- 
vide neutron radiography with an unparalleled advantage over 
other radiographic techniques, such that neutron radiography 
is currently used in the diagnosis of a wide variety of indus- 
trial materials and is of great importance in the metal and 
aerospace industries [1-3]. Simply put, neutron radiography 
takes advantage of the transmission of neutrons and the ab- 
sorption of neutrons by the sample being evaluated, among 
other properties. By comparing the imaging of the detector 
in the presence or absence of the sample, the physical con- 
dition of the sample can be determined. However, there are 
inherent errors in this technique, mainly due to the scatter- 
ing of neutrons by the sample, absorption of neutrons by the 
detector material, and attenuation of neutrons during trans- 
mission. In response to this problem, two main correction 
methods have been proposed: the first is to modify the en- 
ergy spectrum of the neutron source to produce more realis- 
tic final detection results [4, 5], and the second is to intro- 
duce various functions to correct for the effect of scattering 
in neutron radiography [6, 7]. Energy spectrum correction 
is proposed mainly to counter the neutron energy spectrum 
hardening problem. This problem becomes more significant 
as the thickness of the sample to be measured increases to 
a point at which energy spectrum correction becomes neces- 
sary. However, in a broader context, the scattering effect is 
more direct and important than the energy spectrum harden- 
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ing effect. Therefore, this paper discusses energy spectrum 
correction as a secondary correction after scattering correc- 
tion. For example, the 14 MeV fast neutron radiography 
(FNR) technique [20], which has been extensively studied in 
recent years, has promising applications with high contrast 
sensitivity and other features. However, it has difficulty deal- 
ing with the impact of scattering effects, which makes this 
technique much less effective and restricts its use in prac- 
tical applications. There are two popular methods for en- 
ergy spectrum correction: the maximum-likelihood expec- 
tation—maximization (MLEM) method, and the quantitative 
neutron photographic analysis (QNPA) method. The MLEM 
method is mainly based on least squares and maximum- 
likelihood estimations to perform statistical correction [21], 
whereas the QNPA method uses the characteristics of neutron 
decay to perform iterative correction of the neutron distribu- 
tion in order to correct the energy spectrum. Therefore, this 
paper discusses energy spectrum correction as a secondary 
correction after scattering correction, with the optimization 
of the scattering effect as the primary objective. For the scat- 
tering correction, the "point scattered function" has been pro- 
posed [8]. It is essentially the probabilistic density distribu- 
tion function of the scattered neutrons detected by the de- 
tector after the neutron beam passes through the sample to 
be measured. However, it is not easy to apply directly be- 
cause of its complex form and the many parameters it con- 
tains. Therefore, scattering correction models based on a va- 
riety of approximate point-scattered functions have been pro- 
posed [9, 10]. In recent years, materials science has advanced 
rapidly, and an increasing number of neutron detectors can 
use very thin detection layers as a converter to optimize neu- 
tron radiography [11-13]. Thus, we propose a new approx- 
imation of the point scattered function called the "thin-film 
approximate point scattered function". This new model has 
a significant computational advantage and the error is suffi- 
ciently low for practical applications. In the following, we 
derive the thin-film approximate point scattered function from 
the conventional point scattered function and then perform 
numerical simulations using software such as Monte Carlo 
N-particle transport (MCNP) and Mathematica to evaluate the 


accuracy of the model and to perform error analysis. 


I. THIN-FILM APPROXIMATE SCATTERING MODEL 
A. Thin-film approximate point scattered function 


We first review the basic theory of the point scattered func- 
tion, which is based on four basic assumptions: 

1. The scattered neutrons are homogeneous. 

2. The neutrons are scattered only once in the sample. 

3. The source and scattered neutrons have the same atten- 
uation coefficient in the sample. 

4. The sample width is large enough to ensure that the scat- 
tered neutrons penetrate the sample and reach the detector. 

These are reasonable assumptions. Because neutrons are 
electrically neutral, each should be homogeneous after scat- 
tering. The vast majority of neutrons also scatter only once 
in the sample and therefore the need for multiple scattering 
corrections is not significant. 

Fig. 1 shows a schematic of neutron radiography and the 
angular response of the detector. 
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Fig. 1. Schematic diagram of neutron radiography and angular re- 
sponse of the detector 


In Fig. 1, r is the coordinate parameter of the detector, d 
is the thickness of the detector converter, 0 is the scattering 
angle, D is the distance from the sample to the detector, T is 
the thickness of the sample, t is the location where scattering 
occurs in the sample, and R and / are both geometric parame- 
ters. The following relationship exists between the geometric 
parameters. 


D+t 
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Then, for the detector coordinate r, the incident neutron 
flux ®sc (E) and the scattered neutron flux of the detector 
®,, (E) can be expressed in the following form. 
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Here, ®o (E) is the neutron flux of the source, X, (E) and 
E; (E) are the scattering macroscopic scattering cross section 
and the absorbing macroscopic scattering cross section of the 
sample, respectively. “ig; (E) is the absorption cross section 
of the detector material. 


We have to take into account the neutron attenuation effect, 
where a single neutron is detected with a probability of 


PHTt= eo Ud/ cos 0 (3) 


The neutron flux at the detector surface is 
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Then, according to the definition of the point scattered func- 
tion we have 
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This function is quite complicated and inconvenient to use. 


(5) 


Many studies related to scattering correction have been con- 


ducted to approximate this function. 


For a very thin detector, it can be considered that d — 0. 
Because the series expansion of the exponential function 
e—~¢ = 1 — Yd + O [2] involves a second-order infinitesimal 
term O [2], we can use this approximation by considering that 
the main and first-order terms occupy 99% of the exponential 
function (i.e., d satisfies 1 — Nd = 0.99e~~%). This implies 
that the following approximation can be applied. 
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Using the geometric relations and the adjustment order, and 
writing the new function as TPScF (r), we obtain 
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The integral in the numerator of Eq. 7 with respect to ¢ can 


be written as 
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This integral can be calculated numerically using the Newton- 
Cotes method after the associated parameters are given. Fi- 
nally, the thin-film approximation of the model is 
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Here, the integral of F is related to the specific energy spec- 
trum. 


B. Pre-processing of the model 


To better apply this model, we need to pre-process the 
model, as described in this section. We focus on Eq. 8 and ap- 
ply the Newton-Cotes method to obtain a numerical solution 
[14] using 


F(E,r) ~ (E+ED sE on— syn} 42 y rie.ani), (10) 


where h represents the step size and m is obtained by solving 
the following equation. 


(2n — 1) hm =T,m = [hn] (11) 


Here, [x]represents the Greatest Integer Function, which re- 
turns the value of the greatest integer, which is less than or 
equal to x. The error that results from Eq. 10 is given by the 
following equation. 
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Eq. 10 can be calculated very quickly once the determined 
thickness T of the sample and the absorption cross section 
X; (E) are known. 

In addition, the integration of TPScF (r) is difficult. 
Therefore, the approach we adopt is to divide the integration 
interval into many small intervals as follows. 
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Here, E € (Eo, Em) is the range of values of the energy 
spectrum of the neutron source. 


n=1 


Thus, we consider the model from a numerical perspec- 
tive, making it less abstract and more computable in concrete 
terms. Next, we use examples to demonstrate the validity of 
the model and then perform an error analysis. 


II. APPLICATION EXAMPLE 


A. Numerical simulation of scattering correction for water 
and polyethylene 


We construct a simple neutron radiography experiment. 
The neutron converter is made of 100jm-thick °Li material, 
and cylindrical metallic aluminum with a radius of 0.6cm 
and a height of 1.1cm is used as the sample, which is lo- 
cated 0.8cm away from the detector. The experimental setup 
is shown in Fig. 2. 

As the neutron source in the experiment, we use thermal 
neutrons ranging in energy from 0.0leV to leV and whose 
energy spectrum satisfies Maxwell’s spectrum. The neutron 
distribution on the converter can be obtained by importing the 
relevant parameters into MCNP. The results are represented 
as density maps for PScF (r) and TPScF (r), in Fig. 1. 

In the figure, r indicates the distance from the detection 
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Fig. 2. Neutron radiography experiment setup. 
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Fig. 3. Results of (a) PScF (r) and (b) TPScF (r) for the simulation 
experiment 


point to the center of the detector plane. Fig. 3(a) repre- 
sents the simulation results using PScF (r) as a point scat- 
tered function, and Fig. 3(b) represents TPScF (r) for com- 
parison purposes. The red areas in Fig. 3 indicate heavy neu- 
tron scattering, whereas the blue areas indicate very light neu- 
tron scattering. The boundary of the red area lies roughly at 
0.6cm, which is generally in line with the actual condition of 
the sample. Note that the boundary between the two areas is 
not particularly clear. The causes behind this are the ideal- 
ization of the neutron energy spectrum and the haphazard na- 
ture of the scattering in the sample. However, there are well- 
established methods to solve this problem, mainly by correct- 
ing the energy spectrum and by using computer algorithms to 
make the edges clearer, thus making neutron imaging more 
representative of the actual experimental conditions [15, 16]. 
Comparing Fig. 3(a) with (b), it can be seen that TPScF (r) 
gives an overall smaller result than PScF (r), which indicates 
that an overall error exists between them. The Relative Error 
Function is therefore defined as 


|TPScF (r) — PScF (r)| 
= PScF (r) 


(14) 


For the experiment described in Fig. 2, the image of the 
relative error function is mapped in Fig. 4 


This result indicates that the error is very small in the cen- 
tral area of the screen and is larger along the edge areas. 
The average error over the entire detection screen is approx- 
imately 0.306. However, this model is valid because of the 
similarity between Figs. 3(a) and (b), so, we further review 
the validity of the model through error analysis below. 
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Fig. 4. Results of the relative error distribution between PScF (r) 


and TPScF (r). 


B. Error analysis 


The error function presented in Fig. 5 resembles an expo- 
nential function in the radial direction, so we can introduce an 
exponential function to adjust for the error. First, we define 


NPScF (r) = TPScF (r) + Ce7*” (15) 


as the new approximate point scattered function. By fitting 
the relative error data, the values of the parameters C’, a can 
be obtained, and then we have 


NPScF (r) = TPScF (r) + 8.7646 x 107~°e° 398°", (16) 


Then, to study the error relative to the original scattered func- 
tion, the distribution is plotted in Fig. 5. 


Fig. 5. 
NPScF (r). 


Relative error distribution between PScF (r) and 


The figure indicates that the relative errors are at a very low 
level in most regions, and a maximum value of 17% and a 
mean value of 2.5% exist across the entire region. This result 
indicates that a further improvement can be made to the ap- 
proximation model, which validates that NPScF (r) is a very 
effective point scattered function for the thin-film approxima- 
tion. 


IV. NEUTRON IMAGING CORRECTION 
A. Quantitative neutron photographic analysis method 


Owing to the nature of the interaction between neutrons 
and matter, the results of neutron imaging are often very in- 
accurate in actuality. We can reduce the error by using more 
advanced experimental equipment, and/or we can improve the 
imaging technique using algorithms to correct an image. One 
such method, the Quantitative Neutron Photographic Analy- 
sis Method uses algorithms to improve imaging technology. 

Quantitative neutron photographic analysis uses the expo- 
nential neutron decay pattern, 


Teample (i; j) = Io (i, j) e FATE), (17) 


where J is the intensity of the neutron, (i, j) is the prime point 
in the image, and X 4 is the thickness of the sample. The 
distribution of substances can be determined using Eq. 17. 
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This is the principle behind quantitative neutron photographic 
analysis method. 

Numerous experimental studies have shown that quantita- 
tive neutron photographic analysis is influenced by three main 
factors. 

(1). Collimation Ratio 

The collimation ratio is the ratio of the distance from the 
neutron source to the sample with respect to the diameter of 
the output aperture. Because neutrons are electrically neutral, 
they cannot be made parallel by focusing the neutron beam 
through the electromagnetic field, and therefore it is not pos- 
sible to achieve full collimation of the neutron beam using a 
collimator. 

(2). Scattered Neutrons 

The presence of substances with high scattering cross sec- 
tion (e.g., water and polyethylene) in the sample can cause a 
large number of scattered neutrons to reach the detector, thus 
distorting and blurring the neutron image. The solution is 
generally to increase the distance between the sample and the 
detector. 

(3). Neutron Energy Spectrum 

Because incident neutrons are not mono-energetic and neu- 
trons having different energies interact with matter having 
different cross sections, the images and results obtained us- 
ing a single cross section parameter can produce large errors. 
In addition, the energy spectrum of the neutrons scattered by 
the sample may have changed from that of the neutron source 
owing to attenuation during propagation. 

The focus of this study is the correction of neutron images 
by applying a correction for scattered neutrons. To further 
investigate the optimization of the proposed model with re- 
spect to neutron imaging, the current iterative algorithm for 
the quantitative analysis of scattered samples is presented in 
Fig. 7. 


In the algorithm, Isu is the normalized neutron intensity, 
Isc is the scattered neutron distribution function, and k is the 
number of iterations. The criterion equation 
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is used to select a convergence coefficient £ that terminates 
the iterative calculations. The upper criterion actually repre- 
sents the relative change in the transmitted neutron intensity 
calculated for two adjacent iterations. As the algorithm con- 
verges, the value of the change should become increasingly 
smaller, gradually converging toward zero. This algorithm 
has been tested in a large number of neutron scattering exper- 
iments and is highly reliable. This iterative algorithm can be 
applied to the thin-film point scattering function model devel- 
oped above. 


<e (19) 


B. Optimization algorithm based on thin-film approximation 
models 


It should be noted that in the iterative algorithm (Fig. 7), 
PScF (r) is calculated based on the sample thickness, which 
is later brought into the calculation of the neutron inten- 
sity prior to the iterative calculations. However, if we use 
NPScF (r) instead of PScF (r), the iterative convergence can 
be greatly accelerated, as NPScF (r) is equivalent to an ap- 
proximate numerical integration of PScF (r). When we ap- 
plied the iterative algorithm to the experiments in Section 
III.A; using this substitution only three iterations were re- 
quired to obtain satisfactory results, as indicated in Fig. 6. 

The results obtained show the reliability and validity of the 
iterative algorithm on the one hand, and the faster correction 
of imaging results using NPScF (r) instead of PScF (r) on 
the other. The iterative calculation using PScF (r) requires 
five iterations to converge to something closer to actuality 
and takes approximately three times as long to run as the cal- 
culation using NPScF (r). Comparing result (a) using the 
NPScF (r) function with result (b) using the conventional 
PScF (r), it can be observed that (b) produces poor imaging 
results, is slower to converge, and exhibits greater distortion 
at the edges, whereas (a) produces good imaging results and 
slight distortion at the edges. 

We followed up with similar simulations for trapezoidal 
objects and cones and found that using NPScF (r) instead 
of PScF (r) for iterative calculations not only improved the 
simulation, but also improved the speed of the algorithm. 

Moreover, in recent years, a “3D detector” has been intro- 
duced for use in neutron photographic techniques [18, 19]. It 
adds 3D microstructures such as trenches or holes to a pla- 
nar detector and fills them with a neutron converter, which 
can significantly improve the detection efficiency. For this 
class of detectors, the model proposed in this paper is also 
applicable, and together with the algorithm in Fig. 7, yields 
optimized results. However, the relationship between the ge- 
ometric parameters of the planar detector (Fig. 1) must be 
recalculated. That is, to apply the proposed model to the 3D 
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Fig. 6. Scattering imaging results using ((a) NPScF (r), and (b) 
PScF (r)) 


detector, the relationship between geometric parameters must 
be solved again. A mathematically feasible approach is to 
use the symmetry of the 3D detector microstructure to solve 
by parts. For example, if the 3D microstructure is axisym- 
metric, then it can be solved by the radius with the center of 
symmetry being the center of a circle, determining the geo- 
metric relationship at a fixed radius, and then solving using 
the method given in Section 2. 


C. Effects of energy spectrum hardening 


For the thermal neutron spectra considered in this study, 
the cross sections interacting with matter approximately obey 


the 1/v law. As neutrons penetrate the sample, low-energy 
neutrons decay more rapidly than high-energy neutrons, ulti- 
mately increasing the proportion of higher-energy neutrons in 
the neutron energy spectrum, a phenomenon known as energy 
spectrum hardening. 

The current solution to this problem is to define an "effec- 
tive neutron cross section to correct for the energy spectrum. 
Considering that the macroscopic cross section is no longer a 
constant, the neutron transmittance is no longer simply expo- 
nential in relation to the thickness of the sample, even in the 
absence of scattered neutrons. Therefore, the energy response 
of the detector should also be considered. We can define an 
effective cross section based on the ratio of the incident neu- 
tron intensity to the transmitted neutron intensity as 


Ísu L e BelT)T, (20) 
Io 
where — Xop (T) represents the effective cross section for a 
sample of thickness T. 
The intensity of the neutrons detected by a detector is typ- 
ically expressed by the following equation: 


R= | ð (E) (1 — g mial) dE, 21) 
E 
and the transmitted neutron intensity is 


TE Í D (E) eT (1 -e> d) AE. 02) 
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According to the 1/v law, which is approximately satisfied by 
thermal neutrons, the conventional neutron cross section can 
be expressed as 


(23) 


where Ei, œ 0.025eV is the energy of a thermal neutron. 
Considering the detector energy response as a whole, £ (E) = 
1 — e 7 Mact(E)d ~ Yay (E) - d, and combining the above 
equations yields 
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(24) 
However, owing to the difficulties of performing integration 
operations and in obtaining relevant data, use of a simplified 
approximation of the above equation is desirable. The ap- 
proach taken in [17] is to apply a linear fit, but we believe 
that the error in doing so would be very large. Instead, we 
can adopt a numerical integration approach, as described in 
Section 2.2 to approximate the above equation. Again apply- 
ing the Newton-Cotes method, let 
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Fig. 7. Neutron radiography quantitative analysis iterative algorithm flow chart 
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where h is the set step size, and note that ®(£) = 0 at both 
E = 0 and E = oo. With the introduction of the effective 
cross section, the use of an effective cross section to correct 
the energy spectrum can make the decomposition lines of the 
spectrum more visible, so that the edges of the sample or large 
variations in density are more clearly reflected in the image. 

However, it would be preferable to make certain modifi- 
cations to NPScF (r), i.e., Eq. 15, to achieve an improved 
energy spectrum correction. 

Maxwell’s spectrum can be formally defined as ® (E) = 
aE exp (—E/b), and the hardening of the energy spectrum is 
essentially an excess retention of fast neutrons and an excess 
depletion of slow neutrons, and is therefore reflected in an 
increase in the parameter b. Considering the specific form of 
Maxwell’s spectrum, an increase in parameter b is equivalent 
to an increase in the background temperature. Thus, we can 
perform the following transformations to account for these 


factors. 
2mm 2 >E 
(E) — ®) (E) = kT” ry? 4/ a R(T+X) (26) 
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The temperature-corrected energy spectrum is then used in 
the step calculations. Parameter A can then be determined by 
fitting the specific energy spectrum to Eq. 26. 

Substituting the corrected energy spectrum into Eq. 13 and 
replacing the mean scattering cross section with the effective 
scattering cross section in Eq. 25 yields 


TPScF (r) 
_ Eii a (Ei) Eae (Ei) Det (Ei) F (Bar) 27) 
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Subsequently, we can correct for the error by applying error 
analysis, as discussed in Section III.B, by choosing a suitable 


exponential function based on the symmetry of the spatial dis- 
tribution. 


V. CONCLUSION 


In this study, an approximate point scattered function is in- 
vestigated for the correction of neutron photographic imag- 


ing. The modeling is carried out in Section II.A, whereas 
Section II.B explains how the algorithm is used to carry out 
the calculations. 

In the simulation experiment described in Section III, the 
results shown in Fig. 3 demonstrate the validity of the model. 
The model was then corrected for the relative error results 
shown in Fig. 4. This correction is reflected in the addition 
of an exponential function to adjust for the error (see Eq. 15), 
and the parameters are determined using the fitting method 
to finally obtain Eq.16 to produce satisfactory simulation re- 
sults. Fig. 5 shows the validity of the correction, indicating a 
mean relative error of 2.5% for the entire neutron converter. 

In Section IV, we apply the model to the popular quanti- 
tative neutron photographic analysis algorithm and find that 


[1] Brenner C. M., Mirfayzi S R, Rusby D R, et al. Laser-driven x- 
ray and neutron source development for industrial applications 
of plasma accelerators. Plasma Phys. Contr. F. 58(1): 014039 
(2015). doi: 10.1088/0741-3335/58/1/014039 
Ersundu A E, Büyükyıldız M, Ersundu M Ç, et al. The 
heavy metal oxide glasses within the WO3-MoO3-TeO2 sys- 
tem to investigate the shielding properties of radiation ap- 
plications. Prog. Nucl. Energ. 104: 280-287 (2018). doi: 
10.1016/j.pnucene.2017.10.008 
[3] Wang B, Zhong S, Lee T L, et al. Non-destructive testing 
and evaluation of composite materials/structures: A state-of- 
the-art review. Adv. Mech. Eng. 12(4): 1687 (2020). doi: 
10.1177/16878 14020913761 
[4] Liu S Q, Biicherl T, Li H, et al. Corrections on energy spec- 
trum and scatterings for fast neutron radiography at NEC- 
TAR facility. Chinese Phys. C, 37(11): 118201 (2013). doi: 
10.1088/1674-1137/37/11/118201 
Hachouf N, Kharfi F, Boucenna A. Characterization and 
MCNP simulation of neutron energy spectrum shift 
after transmission through strong absorbing materi- 
als and its impact on tomography reconstructed image. 
Appl. Radiat. Isotopes 70(10): 2355-2361 (2012). doi: 
10.1016/j.apradiso.2012.06.017 
Hassanein R, Lehmann E, Vontobel P. Methods of scat- 
tering corrections for quantitative neutron radiography. 
Nucl. Instrum. Meth. A 542(1-3): 353-360 (2005). doi: 
10.1016/j.nima.2005.01.161 
Raventós M, Lehmann E H, Boin M, et al. A Monte Carlo ap- 
proach for scattering correction towards quantitative neutron 
imaging of polycrystals. J. Appl. Crystallogr. 51(2): 386-394 
(2018). doi: 10.1107/s1600576718001607 
Kardjilov N, De Beer F, Hassanein R, et al. Scattering cor- 
rections in neutron radiography using point scattered func- 
tions. Nucl. Instrum. Meth. A 542(1-3): 336-341 (2005). doi: 
10.1016/j.nima.2005.01.159 
[9] Hassan M H. Point Scattered Function (PScF) for fast neu- 
tron radiography. Nucl. Instrum. Meth. B 267(15): 2545-2549 
(2009). doi: 10.1016/j.nimb.2009.05.062 
[10] Xu X, Xie X, He H, et al. Imaging objects through scattering 
layers and around corners by retrieval of the scattered point 


[2 


[e 


[5 


faa 


[6 


= 


[7 


uo 


[8 


fiat 


replacing the conventional point scattered function with the 
proposed approximate point scattered function results in bet- 
ter neutron imaging, faster convergence of the algorithm, and 
a significant reduction in the number of iterations. In ad- 
dition, we investigated the effect of energy spectrum hard- 
ening on neutron photographic imaging, concluding that en- 
ergy spectrum hardening is equivalent to increasing the back- 
ground temperature in Maxwell’s spectrum, and therefore a 
temperature correction parameter can be introduced (Eq. 26). 
After replacing the mean scattering cross section with the ef- 
fective scattering cross section and including the temperature 
correction, a new point scattering function is given by Eq. 27. 
Finally, the error analysis method applicable to the proposed 
model is discussed. 


spread function. Opt. Express 25(26): 32829-32840 (2017). 
doi: 10.1364/o0e.25.032829 

[11] FernandezIzquierdo L, ReyesBanda M G, Mathew X, et al. 
Cesium Lead Bromide (CsPbBr3) ThinFilmBased SolidState 
Neutron Detector Developed by a SolutionFree Sublimation 
Process. Adv. Mater. Technol-US 5(12): 2000534 (2020). doi: 
10.1002/admt.202000534 

[12] Tchouaso M T. Efficiency limitation of thin film coated planar 
semiconductor thermal neutron detector. Appl. Radiat. Isotopes 
173: 109716 (2021). doi: 10.1016/j.apradiso.2021.109716 

[13] Avila-Avendano C, Mejia I, Garcia-Lozano R, et al. Inte- 
grated thin-film radiation detectors and in-pixel amplifica- 
tion. IEEE T. Electron Dev. 65(9): 3809-3815 (2018). doi: 
10.1109/ted.2018.2859778 

[14] Sauer T.. Numerical Analysis. Pearson 231 (2006). 

[15] Hassanein, René Karim. Correction methods for the quanti- 
tative evaluation of thermal neutron tomography. ETH Zurich 
(2006). 

[16] Akkaya, Gizem. "The calculation of self-shielding correc- 

tion factors for large samples in 241Am-Be isotopic neutron 

source." Applied Radiation and Isotopes 179: 109990 (2022). 

doi: 10.1016/j.apradiso.2021.109990 

Capelli S. C., Romanelli G.. An effective hydrogen scattering 

cross section for time-of-flight neutron experiments with sim- 

ple organic molecules. J. Appl. Crystallogr. 52(5): 1233-1237 

(2019). doi: 10.1107/s1600576719011592 

[18] Uher J., Fröjdh C., Jakůbek J., et al. Characterization 
of 3D thermal neutron semiconductor detectors. Nucl. In- 
strum. Meth. A: Accelerators, Spectrometers, Detectors 
and Associated Equipment, 576(1): 32-37 (2007). doi: 
10.1016/j.nima.2007.01.115 

[19] Douglas S., Steven L., Bellinger J., et al. Present status of 
microstructured semiconductor neutron detectors, J. of Cryst. 
Growth, 99(379), (2013). doi: 10.1016/j.jcrysgro.2012.10.061. 

[20] Lu C. B., Bao J., Huang Y., et al. Contrast sensitivity in 14 MeV 
fast neutron radiography. Nucl. Sci. Tech., 28(6): 1-8, (2017). 
doi: 10.1007/s41365-017-0228-5. 

[21] Beni M. S., Krstic D., Nikezic D., et al. Studies on unfold- 
ing energy spectra of neutrons using maximum-likelihood ex- 
pectation—maximization method. Nucl. Sci. Tech, 30(9): 1-10, 
(2019). doi: doi: 10.1007/s41365-019-0662-7. 


[17 


(am 


